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The staggered quark action on anisotropic lattices is studied. We carry out numerical 
simulations in the quenched approximation at three values of lattice spacing (a^ ^ = 1 — 2 
GeV) with the anisotropy ^ = Oct/cIt = 4, where and are the spatial and temporal 
lattice spacings, respectively. The bare anisotropy 7f in the quark action is numerically 
tuned through the ratio of meson masses in the fine and coarse directions, and through the 
dispersion relation of a meson, so that the renormalized fermionic anisotropy coincides with 
that of the gauge field. The discrepancy between these two calibration schemes provides 
an estimate of the finite lattice artifact, which is found to be sizable in the range of cutoff' 
explored in this work. We also compute the meson masses using correlators with the wall 
source at the tuned anisotropy parameter. The flavor symmetry breaking effect smoothly 
decreases as /3 increases. The effect of uncertainty in 7f on the meson masses are examined. 
We also discuss a perspective on dynamical simulations. 

§1. Introduction 



O . In numerical studies of lattice QCD, one frequently encounters a case which re- 

quires a fine lattice spacing in the temporal direction while does not comparatively in 
the spatial directions. Increasing the lattice cutoffs in all the four directions severely 
^ , increases the computational cost, since it increases at least in proportion to the vol- 

^ I ume of the lattice, and in fact more rapidly in particular in dynamical simulations.^^ 

A solution is provided by anisotropic lattices, on which the temporal lattice spacing 
is finer than the spatial ones. The technique is useful in various fields of the lattice 
QCD simulation: At finite temperature, a large number of the degrees of freedom in 
the Euclidean time direction leads a large number of Matsubara frequencies, which 
is efficient for calculations of the equation of state'^^ and for analyses of temporal 
correlation functions of hadrons.^-*"^^ The large temporal cutoff is important for rel- 
ativistic formulations of heavy quark on the lattice. ^^"^^^ It is also efficient when the 
signal-to-noise ratio deteriorate quickly, as in the cases of glueballs,^^' negative 
parity baryons,^^^ and the pion scattering length. 

On the other hand, on anisotropic lattices one has additional parameters in 
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the actions which control the anisotropics of the fields. In general these anisotropy 
parameters should be tuned numerically. Inappropriate tuning of the parameters 
breaks the rotational symmetry of the lattice, and may lead unphysical results. The 
uncertainties in the tuning of the parameters bring additional errors into observed 
quantities. For precise calculations, one needs to tune anisotropy parameters with 
good statistical accuracy and to control the systematic errors in the continuum ex- 
trapolation. In this paper, we investigate this calibration process in detail for the 
staggered quark action. 

To explain the situation we are faced with in more detail, let us focus on a study 
of hadrons at finite temperature, which is a main motivation of present work. To 
investigate the hadron structure at T > 0, one needs to treat the hadron correlators 
in the temporal direction.^)' ^^'^''^ Because the lattice size in the temporal direction 
is limited to NfUr = l/T, where T is the temperature, near and above the critical 
temperature the number of degrees of freedom is severely limited without introducing 
the anisotropic lattice. For example, the number of degrees of freedom is significant 
for reliable extraction of the spectral function from the lattice data.^-*'^)'^^) 

So fax the studies of hadrons at finite temperature have been performed mainly 
on quenched lattices. Dynamical simulations are manifestly important, as exhibited 
by the fact that the order of phase transition changes as the number of dynamical 
quark flavors varies. As the quark action wc adopt the staggered action, which 
has several advantages over the Wilson-type quarks:^^^ Firstly, the staggered action 
retains the remnant of the chiral symmetry. This is important to investigate the role 
of the chiral symmetry near the phase transition. Secondly, one can explore smaller 
quark mass region than the cases with the Wilson-type formulations. Thirdly, the 
cost of numerical simulation is much more economical than the other formulations. 
A disadvantage lies in its complication in the flavor structure. To circumvent the 
effects of the flavor symmetry breaking, improved versions of the staggered fermion 
have been developed. 

In this work, however, we adopt the simplest version without any improvement. 
This is because an improvement adds the anisotropy parameters which is to be tuned 
in general nonperturbatively. At this first stage of development of the anisotropic 
staggered action, we concentrate on the most significant effect of the anisotropy on 
the spectrum. For the gauge field we adopt the standard Wilson plaquette action, 
which has the discretization errors in the same order as the quark action we adopt. 
The same combination of actions have also been used in Ref.^''-* for a dynamical 
simulation, while their discussion on the systematic errors and statistical precision 
were not sufficient for our present purposes. 

The dynamical anisotropic lattices arc more involved than the quenched case, 
since one needs to tune the anisotropy parameters for the gauge and quark fields 
simultaneously. So far there is no systematic investigation of anisotropic staggered 
quarks even in the quenched approximation. In this paper we therefore concentrate 
on the quenched approximation and investigate the properties of the staggered quark 
on anisotropic lattices. 

The numerical simulations are performed on quenched lattices with three lattice 
spacings, at fixed renormalized anisotropy ^ = 4. These scales cover the range of the 
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spatial lattice cutoff = 1 — 2 GeV. We apply two calibration procedures with 
different definitions of fermionic anisotropy S^p: with the ratio of the masses in the 
fine and coarse directions, and through the meson dispersion relation. The former 
has an advantage in statistical precision, while the latter can be used even in the 
heavy quark mass region. Differences between the results of these two procedures 
signal the finite lattice artifacts. For the reason of statistical fluctuations, we use only 
the pseudoscalar channel in the calibration. Other mesonic channels are observed in 
the second part of the paper by employing the wall source, and compared with the 
results on the isotropic lattices. 

This paper is organized as follows. The next section briefly summarizes the 
anisotropic staggered action. Section |21 describes the calibration procedures. The 
numerical result of the calibration is presented in Sec. 0] Section [S] describes the 
spectroscopy with the wall source. Section El is devoted to our conclusions in the 
quenched simulations. In the last section, we give our perspective toward the dy- 
namical simulations with the staggered quarks. 

§2. Staggered quark on anisotropic lattices 

2.1. Anisotropic lattice actions 

The gauge field is described with the standard Wilson gauge action, 

3 1 r n 3 
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ReTr Ui^ix 



, (2-1) 



where /3 = Q/gQ is the bare coupling and 7g the bare gauge anisotropy param- 
eter. The plaquette variable U^y{x) is defined with the link variable U^{x) ~ 

^iga^A^,(x) £SU(3) as 



U^,,{x) = U^{x)U,{x + /l)C/t(x + y)Ul{x). 



(2-2) 



A lattice site is labeled by an integer vector x whose component x^ is in units of 
lattice spacing a^, where ai = 02 = as = Co- and 04 = a,-, /i is a unit vector in the 
/i-th direction. 

The staggered quark action on an anisotropic lattice is defined as 



Sf = ^x{x)K{x,y)x{y), 



(2-3) 



x,y 



K{x, y) = 5x,y - /tfj ^ mix) Ui{x)6, 



x+i,y 



U}(x-i)5 



1=1 



lFt^am{x) Ui{x)5^_^i -Ul{x-A)5. 



z— 4,3/ 



(2-4) 



where is the bare anisotropy of the quark field, = l/2mg with ruq the 
bare quark mass in spatial lattice units. rj^{x) is the staggered phase, = 
(_l)^i+-+^M-i, with r/i(x) = 1. 
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The staggered quark fields x a-nd x have only the color components. The 4-spinor 
field with four degenerate flavors, il)a{X), where X labels the space-time position in 
units of 2a^, are represented as a linear combination of x{'^) as^^^'^^^ 



W = 7f E(^pW^p(^)' (2-5) 



rp = 7f72'^7f^7f, (2-6) 
Xp{X) = x{2X + p), (2-7) 

where p labels a position in the hypercube and its components are or 1. Therefore 
the quantum number of a quark bilinear is specified by its spin-flavor structure. 

2.2. Meson correlators 

The meson correlators are composed of the quark propagator, 

S{x, y) = {x{x)x{y)) = K-\x, y), (2-8) 
which is obtained by solving a linear equation 

^K(z,x)5(:e,0) =6(z)5,,,0, (2-9) 

X 

where the vector h{z) is a source field. (The color index is omitted.) As h{z), we 
select the following two cases; (i) the point source, bp{y) = 6^^^, and (ii) the wall 
sources, beiy) = 1 and bo{y) = (— l)2^i+?'2+y3_ jj^ the latter case, the gauge must be 
fixed. 

In this paper, we treat single time-slice correlators, namely the sink operators at 
each t are defined only in a spatial plane on a single time slice. Such a correlator 
contains two modes, one monotonously and the other oscillatingly decay. The latter 
has the opposite parity to the former. In the following, the correlators are labeled 
by their non-oscillating channels. 

The meson correlators composed of the quark propagator with the point source 
are used in the calibration (Sec. EJ)- Since we do not fix the gauge at this stage, 
the correlators must be constructed with the gauge invariant ingredients. The pseu- 
doscalar and vector correlators are represented as 

Mpsit) = Y,\Spix)\^ (2-10) 

X 

Mv,{t) = T.i-r^\Spix)\^ U = 1,2,3), (2-11) 

X 

where Sp{x) is the quark propagator with the point source at the origin. These 
correlators have spin-fiavor structures 75 75 and jj 7^ , respectively. While these 
representations are for the vanishing meson momentum, the momentum insertion is 
straightforward. In the case of vector correlator, the finite momentum state needs 
some care, since the operator mixes with the fourth component. However we do not 
discuss this problem further, since we decide to use only the pseudoscalar correlator 
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in the calibration because of large fluctuation in the vector channel. We also measure 
the correlators in the z-direction to define the anisotropy through the ratio of masses 
in the fine and coarse directions. Rewriting Eqs. ()2-lU|) and (|2-lll) for the correlators 
in z-direction is also straightforward. 

The correlators composed of the quark propagators with wall sources are as 

(2-12) 
(2-13) 
(2-14) 
(2-15) 
(2-16) 
(2-17) 

The spin- flavor structures of these correlators are; vr: 75 75, vf: 7475 (gi 7475, ir^: 
75 ® 7573, 71-3: 7475 (gi 747573, pi- 73 72, and p^: 7374 7274. In addition to 
these quantum numbers, each correlator contains the opposite parity channel as the 
oscillating mode. The correlator (|2-12() has the same quantum number as Eq. ()2-10() . 
On the other hand, the quantum number of the correlators H2-16|) and H2-17() are 
different from that of Eq. H2-ll() at finite lattice spacing. 
At large t, the correlator approaches a form^*^) 

M{t) Zexp(-mt) + (-)*Zexp(-mi). (2-18) 

To extract the meson masses, we fit the numerical data of the meson correlator with 
the wall source to this form after adding the contribution via the temporal boundary. 
In the calibration, only the first term in Eq. is retained since the contribution 

from the second term quickly disappears in the pseudoscalar channel. 

§3. Calibration procedures 

3.1. Calibration schemes 

The anisotropy parameters 7g and should be tuned so that the physical 
isotropy condition, Ccilc'^F) ~ ^Fi^c^^p) ~ ^' ^olds at each set of (3 and rriq, 
where icinGilF) and ^FildF) are the renormalized anisotropics defined through 
the gauge and fermionic observables, respectively. In quenched simulations, one can 
firstly determine 7^ independently of 7ir, and then tune 7^ for fixed 7^. In contrast, 
7g and need to be tuned simultaneously in dynamical simulations. 



follows. 

M^t) = [sl{x)Se{x) + Slix)Soia 
M^{t) = ^(-)-i+-2+-3 [sUx)So{x) + Sl{x)Se 

X 

MnAt) = [^ei^ + - Slix + 3)So{x) 

X 

Mi,{t) = ^(-)-i+-2+-3 ^slix + 3)So{x) - Slix + 3)5e 

X 

MpA{t) = [slix + 2 + 3)Seix) + slix + 2 + 3)Soix) 

X 

^Pi^*^ = Yi-T'^""'^"' Hi^ + 2 + 3)5o(x) + slix + 2 + 3)5e(a 
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In this work, we use the values of the gauge parameters [5 and 7^ whose renor- 
mahzed anisotropy, ^g, has been determined in good accuracy. The values of the 
parameters used will be described in detail in Sec. 14.11 For the quark field, we em- 
ploy the following two calibration schemes accordingly to the definitions of fermionic 
anisotropy £,f- 

• Mass ratio scheme: In this scheme, is defined as the ratio of meson masses 
in the temporal and spatial directions, 



(^) / (i) 



(3-1) 



We use only the pseudoscalar channel as already noted in the previous section 
for the statistical reason. 

Dispersion relation scheme: An alternative definition of (,f makes use of the 
meson dispersion relation,^) 



(3-2) 



E(p) and mu are in temporal lattice units and p is in spatial lattice units, and 
hence appears. On a finite volume lattice, pi = 27rnj/Lj {i = 1,2,3), where 
Li is the lattice size in the i-th direction. Eq. H3-2|) requires that the rest mass 
equals the kinetic mass. 

3.2. Free quark case 

The mass dependences of ^^"^ and ^j^^^ for the free quark propagator give us a 
guide to analyze the results of the numerical simulation. The free quark propagator 
of the staggered quark field 1)2 -7^ in the momentum space is'^^-* 



Spp'ip) 



i Ej ^pp'iP) sin(Pj/2) - ir^p'iphF sin(p4/2) + mqS^ 



pp 



^pp 



where 



sm'{pj/2) + 72 sm'{p4/2) + 

Jp{p~p')/2i 



The propagator in the t-direction, 

has asymptotic behavior Sppi{p,t) oc e^^'*'^^*, with E^'^\p) implicitly given by 

2 



(3-3) 
(3-4) 

(3-5) 



chE^^\p) = 1 + 



■ 2 Pi 2 
sm — + m„ 

2 



(3-6) 



Similarly, for the propagator in z-directions. 



chS(3) (p) = 1 + 2 



i=i,2 



■ 2 Pj , 2 ■ 2 ■P4 . 2 

sm — + 7p sm h m„ 



(3-7) 
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Table I. The lattice parameters. The values of 70 are based on the resuhs of Ref. (/3 = 5.95, 
6.10) and Ref.^°^ (/3 = 5.75). The lattice scale a''^ is set by the hadronic radius ro.^°^'^^^ The 
mean-field values are from Ref.^^-* 'Note that those values at f3 — 5.75 were evaluated at slightly 
different anisotropy, 70 = 3.072. 



p 


7G 




Size 


Nconf 


ro 


a-'{GeV) 






5.75 


3.136 


12^ 


X 24 X 96 


224 


2.786(15) 


1.100(6) 


0.7620(2)* 


0.9871* 


5.95 


3.1586 


16^ 


X 32 X 128 


200 


4.110(23) 


1.623(9) 


0.7917(1) 


0.9891 


6.10 


3.2108 


20^ 


X 40 X 160 


200 


5.140(32) 


2.030(13) 


0.8059(1) 


0.9901 



where p = {pi,p2,P4:)- For a small quark mass, the rest masses. Mi = E{p = 0), are 
represented as 



(4) 



1 



-rUn 



(3) 



IF 
2m„ 



3-m3 + 0(m^^ 



3x 

+ 0{m\ 



Requiring 4*^^ = ikff ^ = 



*{M) 

If 



1 + 
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1 



^2 / rnl + 0{ml). 



(3-8) 



(3-9) 



By differentiating Eq. (|3-6j) by pj, the tuned anisotropy in the dispersion relation 
scheme is obtained as 



*(DR) 

7f 



1 



1 



37 



2^1 + O(m^) 



(3-10) 



Therefore in both the schemes the quark mass dependences of 7^ are of O(m^). As 

the bare quark mass increases, ^*p^'^ increases while j^^^^ decreases. These results 
are useful to understand the results of the numerical simulation. 



4. Numerical result of calibration 



4.1. Lattice setup 

For the quenched anisotropic lattices, several results for the calibration of gauge 
field are available. ^^^^^''^ We adopt three sets of the lattice parameters (/3,7g) at 
the renormalized anisotropy ^ = 4, according to Refs.^^^'^'^^ The parameters are 
summarized in Table These lattices have almost same physical volume. The 
temporal direction is chosen as the fine direction, and the size in the z-direction is 
set being the same physical length as the f-direction. The lattice scales are set by 
the hadronic radius vq.^^^ 

For the larger two values of (3, we adopt the Klassen's result for 7^ which was 
determined within 1% statistical accuracy. These two lattices are almost the same 
as those used in Ref.,^^^ while the size in the z-direction is doubled. 

At /3 = 5.75, we use a recent result of a precision calibration using the hadronic 
radius tq^^^ as a calibration condition.^'''' With the Liischer-Weisz noise reduction 
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Fig. 1. The effective mass plots for — 0.02 and 7f = 2.8 at /3 = 5.75. The left and right panels 
show the effective masses in the t- (fine) and z- (coarse) directions, respectively. 



technique, '^^^ the static quark potential was computed to the level of 0.2% accuracy 
in the fine and coarse directions. Since the value used in this paper was obtained in 
an earlier stage of the work,^*') 7^ is slightly different from the final result quoted in 
Ref.30) 

We also quote the spatial and temporal mean-field values,"^^^ and Ur respec- 
tively, in the Landau gauge determined in Ref.^^^ For (3 = 5.75, the values listed 
in Table ^ are not correct exactly, since they were determined at the different 
(7g = 3.072). This is not a serious problem, since we quote these values just for 
a qualitative comparison of the numerical result of 7^ with a mean-field estimate, 

4.2. Correlator data 

In the calibration, we use the pseudoscalar meson correlators with the point 
source. To perform the calibration along the two schemes described in Sec. |H1 we 
need the following two types of correlators for each set of rUg and 7^: 

• Correlators in the fine (t-) direction at finite momenta. We compute them at 
momenta pi = 27rnj/Lj {i = x,y,z) where n =(0,0,0), (0,0,1), (0,0,2), (1,0,0), 
(1,0,1), (1,1,0), and (1,1,1). Note that = Ly = L,/2. 

• Correlators in the coarse (2-) direction at zero momentum. 

To observe whether the correlator is dominated by a single state and to determine 
a fit range, we observe an effective mass ^i^e// which is defined through 

M(t + 2) ^ exp[-me//(t + 2)] + exp[-me//(A^f - ^ - 2)] ^^^^^ 
M{t) exp[-me//t] + exp(-me//(iVf - t)] 

The definition with lattice spacing 2a is adopted since the correlators are composed 
of monotonously and oscillatingly decaying modes. An effective mass in z-direction 
is defined similarly. Figure ^ displays the effective masses in the t- and z-directions 
at /3 = 5.75, niq = 0.02, and = 2.80. 



9 




0.2 0.4 
(lattice unit) 



0.2 0.4 
(iattice unit) 



Fig. 2. The dispersion relation at /3 = 5.75. The left and right panels display the results at 
(mg,7F) = (0.02,2.8) and (0.10,2.878), respectively. 



Figure ^ shows that around t = 24 and z = 6 the correlators reach plateaus 
where the contribution from the oscillating modes and excited states is sufficiently 
reduced. We choose a fit range for each quark mass, while the same fit range is used 
for all the momenta, and fit the correlators to a single exponential form. 

Taking the ratio of the masses observed in the t- and z-directions, the fermionic 
anisotropy ^^^^ is determined for each input 'jp. 

Figure 121 shows the meson dispersion relation for two quark masses at /? = 5.75. 
E(p)'^ is almost a linear function of for small momenta, and well fitted to a 
quadratic form in the whole measured momentum region. The same tendency is 
observed for all the quark masses and /3's explored. We apply a linear fit to the data 
of E{p)'^ with the smallest five momenta, and determine the fermionic anisotropy 
^^^^ through its slope. As displayed in Fig. El the linear and quadratic fits give 
almost the same results, while the former gives a smaller statistical error. 

4.3. Determination of 

We observe the fermionic anisotropics i^^'' and ^^^"^ for several values of 
at each quark mass. Figure displays typical 7i? dependences of ^i?'s at /3 = 5.75. 
This figure shows that the renormalized anisotropy is well represented as a linear 
function of the bare anisotropy around £,f = ^- At /5 = 5.95 and 6.10, we therefore 
in most cases measure the correlators for two values of 7ir at each quark mass. 

Interpolating ^_f(7f) linearly in 7^? to ^ = 4, we can define a tuned anisotropy 
parameter 7^. To specify the tuning procedure, the 7^ tuned with the mass ra- 
tio (dispersion relation) scheme is denoted by 7^^^^ (7^^^^)- '^^^ result of the 
interpolation is summarized in Table UTI and displayed in Figured 

Figure 131 shows large discrepancies between the results of two calibration schemes 
in the large quark mass region. The differences remain at the massless limit for 
(3 = 5.75 and 5.95, while disappears within errors at /3 = 6.10. This implies that 
the discrepancy is due to the finite lattice spacing artifacts. The analysis of the free 



10 




3.8 



3.6 



Fig. 3. The dependences of £,f's on 7f at /3 = 5.75. The left and right panels are for niq — 0.02 
and 0.10, respectively. The solid lines are the results of quadratic fits just for a guide of eye, 
while in the determination of 7^ the linear fits using the data of near ^ are used. 




Fig. 4. The results of the calibration. The left panel displays the results at /3 = 5.75 and /3 = 5.95. 
The right panel shows the results at /? = 6.10. The data are from Table HTl The solid and dashed 
lines represent the results of fits which are listed in Table UTTl 



quark propagator suggests that the quark mass dependence of 'jp on nig starts with a 
quadratic term near the chiral limit in both the schemes. This tendency is observed 
in the dispersion relation scheme more clearly than the mass ratio scheme. As the 
bare quark mass increases, 7^''^^'' increases while ^'^^^'^ decreases. These quark mass 
dependences are consistent with the results for the free quark case, Eqs. ()3-9() and 
(|3-10p . For ^ = 4, these expressions expect larger quark mass dependence for 7^*'*^'* 
than ^*p^^\ which is also observed in the numerical result. 

The discrepancy between 7^^^^^ and 7^^^^^ at the chiral limit can be explained 
by the O(a^) finite lattice artifacts, such as the flavor symmetry breaking effect, the 
0((ap)^) uncertainty in the dispersion relation scheme, and so on. The latter effect 
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Table II. The results of the calibration. 7^ and 7^ represent the tuned anisotropy param- 
eters in the mass ratio and dispersion relation schemes, respectively. The pseudoscalar quark 
masses at 7J. are determined by interpolation. 



■mq input 7F 7*'"' rn%{-y*p^^^^) 7f"'"'' ^ps(7f'^^ 
P = 5.75 

0.50 3.578, 3.5761, 3.136, 2.80 3.5822(21) 0.46494(17) 2.6777(47) 0.59849(73) 

0.40 3.471, 3.466, 3.136, 2.778 3.4735(23) 0.42055(17) 2.7239(45) 0.51765(61) 

0.30 4.00, 3.353, 3.30, 2.80 3.3590(27) 0.36962(17) 2.7700(47) 0.43311(49) 

0.20 3.219, 3.136, 2.80 3.2308(35) 0.30798(19) 2.8073(60) 0.34376(50) 

0.10 3.136, 3.08, 2.878, 2.80 3.0870(34) 0.22486(16) 2.821(11) 0.2398(10) 

0.05 3.136, 3.00, 2.987, 2.80 2.9891(76) 0.16339(18) 2.8296(95) 0.16979(37) 

0.03 3.136, 3.00, 2.948, 2.80 2.947(10) 0.12867(21) 2.827(10) 0.13254(33) 

0.02 3.136, 3.00, 2.918, 2.80 2.930(13) 0.10630(22) 2.831(10) 0.10845(29) 

0.01 3.136, 3.00, 2.888, 2.80 2.924(28) 0.07620(36) 2.840(12) 0.07765(23) 
/3 = 5.95 

0.50 3.624, 3.578 3.6246(16) 0.45727(13) 2.8513(94) 0.54393(98) 

0.40 3.521, 3.471 3.5198(19) 0.41144(14) 2.897(10) 0.47370(98) 

0.30 3.413, 3.353 3.4113(25) 0.35813(15) 2.945(12) 0.39781(99) 

0.20 3.292, 3.219 3.3012(36) 0.29292(17) 3.000(15) 0.31307(97) 

0.10 3.166, 3.080 3.1958(41) 0.20575(12) 3.032(13) 0.21278(53) 

0.05 3.112, 2.987 3.1460(91) 0.14498(18) 3.043(25) 0.14781(65) 

0.03 3.170, 3.100 3.139(17) 0.11279(28) 3.053(77) 0.1144(13) 

0.02 3.100, 2.918 3.137(20) 0.09258(27) 3.046(80) 0.0940(11) 

0.01 3.150, 3.100 3.131(31) 0.06631(27) 3.05(12) 0.0670(11) 

/3 = 6.10 

0.30 3.21, 3.11 3.4341(26) 0.34436(17) 3.0174(95) 0.38259(96) 

0.20 3.21, 3.11 3.3309(35) 0.27754(16) 3.068(10) 0.29493(74) 

0.10 3.21, 3.11 3.2377(69) 0.18959(19) 3.119(16) 0.19413(62) 

0.05 3.21, 3.16, 3.10 3.199(16) 0.13037(24) 3.158(32) 0.13129(71) 

0.03 3.31, 3.21, 3.16 3.199(29) 0.10009(25) 3.198(46) 0.10010(83) 

0.02 3.31, 3.21, 3.16 3.210(45) 0.08155(29) 3.244(66) 0.08113(95) 

0.01 3.31, 3.21, 3.16 3.250(80) 0.05798(34) 3.261(78) 0.05789(68) 



concerns the assumed form of the meson dispersion relation as well as the fitting 
form for the energies at finite momenta. For the measurement of the mass in the 
z-direction, the meson quantum number is composed of the staggered quark fields 
on asymmetric cubes in a z-plane. This leads to a different manifestation of the 
flavor symmetry breaking effect from that of the mass in t-direction. Although all 
these effects must disappear toward the continuum limit, at finite lattice spacings 
they may cause intricate systematic effects. 

We determine 7^'s at the vanishing quark mass by fitting the data to a linear 
form in or quadratic form in nig, 

1*F = Co + Cimg + C2ml (4-2) 

where Ci is set to zero in the linear fit. The numbers of data used in the fits are 
selected appropriately. The result of the fit is listed in Table llll| and displayed in 
Fig. 13 For both the schemes, the linear form in seems to represent well the data 
in the light quark mass region. Although the quadratic form does not seem to work 



12 



Table III. The result of the fits of the tuned anisotropy parameter to the form l|4-2|l . In the 
case of linear fit, (i is set to zero. In the second column, 'M' and 'DR' represent the mass ratio 
and dispersion relation schemes, respectively. Ndata is the number of data points used in the fit. 







xnass ratio 




disp 


srsion relation 






scheme 


fit 


Ndata 


Co 


Ci 


C2 




5.75 


M 


linear 


4 


2.940(11) 




14.7(11) 


0.092 




M 


quadratic 


4 


2.915(42) 


0.4(19) 


21.(23) 


0.001 




DR 


linear 


5 


2.833(10) 




-1.2(12) 


0.013 




DR 


quadratic 


6 


2.837(13) 


-0.19(29) 


0.2(12) 


0.009 


5.95 


M 


linear 


5 


3.131(14) 




6.4(12) 


0.0007 




M 


quadratic 


5 


3.137(40) 


-0.2(12) 


7.9(78) 


0.0006 




DR 


linear 


5 


3.049(38) 




-1.6(35) 


0.0002 




DR 


quadratic 


6 


3.053(71) 


-0.2(11) 


-0.6(38) 


0.0002 


6.10 


M 


linear 


5 


3.194(28) 




4.3(24) 


0.008 




M 


quadratic 


5 


3.242(81) 


-1.7(21) 


17.(13) 


0.002 




DR 


linear 


5 


3.206(45) 




-8.9(40) 


0.021 




DR 


quadratic 


6 


3.256(69) 


-1.9(10) 


4.6(34) 


0.007 



0.114 p 

0.112 - 

0.110 - 

0.108 - 

0.106 - 



P=5.75 m =0.02 
PS 



2.9 3.0 3.1 

Vf 




Fig. 5. The dependence of mps on 7f at /3 = 5.75. The left and right panels show the results at 
TUq — 0.02 and 0.10, respectively. The solid lines represent the results of quadratic fits. 



properly, the 7^ at the chiral limit is close to the result of the linear fit. 

4.4. Pion mass and the chiral limit 

Here we verify that the pseudoscalar meson mass satisfies the PCAC relation 
with the calibrated anisotropy parameter. Firstly we need to interpolate the pseu- 
doscalar meson mass to 'yp = 1*f- Figure El shows typical dependences of PS 
meson masses at (3 = 5.75. The PS meson mass is well represented by a linear form 
in the vicinity of 7^. The result of the linear interpolation is listed in Table ITTl 

Figure ini displays the niq dependence of the pseudoscalar meson mass squared. 
The proportionality of m?pg to ruq well holds for both the calibration schemes. This 
is also verified by linear fits whose results are listed in Table HVl In either scheme, 
the intercept is quite small and becomes consistent with zero within the error as 
(3 increases. The difference between the two schemes also decreases toward the 
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0.02 



o p=5.75: mass ratio 




• p=5.75: dispersion relation 




□ p=5.95: mass ratio 




■ p=5.95: dispersion relation . ' 




^ P=6. 10: mass ratio ^."'x 




A p=6. 10: dispersion relation m'X 









0.00 0.01 0.02 0.03 

Fig. 6. The dependence of •m^'pg on the bare quark mass m,. The sohd and dashed hnes represent 
the results of linear fits. 



Table IV. The results of linear fits of mpg in at the tuned anisotropy 7J. 



/3 


mass ratio 
mj(0) dm%/dmq 




dispersion relation 
m^(0) drn^/dniq 


x' 


5.75 
5.95 
6.10 


0.000493(38) 
0.000236(36) 
0.000028(43) 


0.5382(17) 
0.4164(22) 
0.3325(16) 


8.14 
0.0407 
0.458 


0.000251(49) 0.5784(36) 
0.00021(20) 0.430(13) 
0.000025(91) 0.3317(70) 


0.535 
0.0214 
0.319 



continuum limit. 

4.5. Summary of calibration 

In the next section, we investigate the meson spectrum of the anisotropic stag- 
gered quark. For this purpose, it is convenient to select one of the results of cali- 
brations as our main result for 7^ and compute the spectrum with this value. The 
differences with the other choices are treated as a systematic uncertainty whose effect 
on the spectrum should be investigated. 

As such a representative, we adopt the 7^ in the massless limit determined by 
the linear fit in in the dispersion relation scheme for the following reasons: As 

apparent in Fig. lU j*p^^^ less depends on the quark mass than 7^^*^^ . In the light 
quark mass region we can use the value of 7^ at the chiral limit instead of the direct 
results on those quark masses. The statistical error of the former can be reduced 
less than those of the latter. 

Estimate of the statistical error of 7^ at massless limit is provided by the sta- 
tistical error of the fit result for ^o- On the other hand, the discrepancy between the 
two calibration schemes gives a typical size of the systematic uncertainty. 

To summarize, we select the following values as a representative result of the 
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calibration: 



6.10 



5.95 



5.75 




2.83(1)(11) 
3.05(4)(11) 
3.21(5)(5), 



(4-3) 



where the first and second parentheses denote the statistical and systematic errors, 
respectively. The latter does not include the quark mass dependence of 7^. At 
P = 6.10, the domination of the statistical error around the chiral limit disable 
us from estimating the systematic error in the same way as other /3's. Prom the 
conservative point of view, we substitute the statistical error for the systematic one, 
since the latter is at most of the same size as the former. 

Finally let us compare the result of calibration with the mean-field estimate. 
The mean-field improvement is performed by replacing the link variables as Ui 
Ui/ua {i = 1,2,3) and — > Ui/uT- The mean-field estimate of 7^ in the chiral 
limit results in 



From the values quoted in Tabled 7^ is determined as 3.0878(8)*^ (/3 = 5.75), 
3.2017(4) (/? = 5.95), and 3.2558(4) {(3 = 6.10). The result of numerical simulation 
approaches to the mean-field estimate from below as (3 increases. This suggests that 
the mean-field value of 7^ provides a good guide in the calibration. This feature is 
helpful in the calibration of anisotropic lattices with dynamical quarks. 



In this section, we compute the meson spectrum using the anisotropic staggered 
quark action tuned in the previous section. As noted at the end of Sec. |11 we 
use the values of quoted in Ea. H4-3() in the simulation and investigate how the 
uncertainties in 7^ affect the spectrum. In practice, we compute the masses at two 
7ir's and linearly interpolate them to 7^ in Ea. (|4-3|) . Simultaneously the response 
of the mass to the change of is obtained. 

The numerical simulation is performed at the same three /3's as in the calibration. 
At /3 = 5.75, the lattice size is the same as in the calibration. At /3 = 5.95 and 6.10, 
we use the lattices with half size in the z-direction, which were used in Ref.'^^^ 
The numbers of configurations are 224, 400, and 200 for (3 = 5.75, 5.95, and 6.10, 
respectively. 

We measure the meson correlators with the wall source. The correlators mea- 
sured were described in Sec. |21 The gauge configurations are fixed to the Coulomb 
gauge. 

The complication in analyzing the spectrum of the staggered mesons arises from 
the oscillating modes contained in the correlators. This makes fit more involved than 
other quark formulations. In extracting the meson masses, we apply the constrained 

*-* At /3 = 5.75 the value has additional systematic error since the mean-field values are obtained 
at slightly different 7^. 




(4-4) 



§5. Spectroscopy on anisotropic lattices 
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Table V. The meson spectra from the correlators with the wall source. The meson masses at 
ruq = are determined by linear fits using the data at lightest three m^'s. 



rria 



p = 0.(0 














u.iu 




U.Z4o i ^ylZ) 


U.Zy4iUt^4U ) 




u.o4ioyt^oo ) 




0.05 


0.16935(12) 


0.17916(32) 


0.22577(43) 


0.22749(69) 


0.28222(94) 


0.2831(16) 


0.03 


0.13211(14) 


0.14245(43) 


0.19110(39) 


0.1941(11) 


0.25546(96) 


0.2557(16) 


0.02 


0.10854(15) 


0.11997(49) 


0.17199(48) 


0.1753(16) 


0.2408(13) 


0.2406(23) 


0.01 


0.07755(15) 


0.09159(60) 


0.14948(68) 


0.1569(22) 


0.2209(22) 


0.2207(34) 


0. 


0.01726(67) 


0.0495(15) 


0.12387(98) 


0.1336(34) 


0.2064(28) 


0.2054(42) 


f3 = 5.95 
0.10 


0.21169( 8) 


0.21571(12) 


0.23084(17) 


0.23136(24) 


0.25809(47) 


0.25816(63) 


0.05 


0.14754(10) 


0.15056(14) 


0.16469(20) 


0.16529(37) 


0.20203(53) 


0.20163(73) 


0.03 


0.11429(13) 


0.11704(17) 


0.13182(24) 


0.13222(50) 


0.17759(74) 


0.1767(11) 


0.02 


0.09380(15) 


0.09660(23) 


0.11250(28) 


0.11296(61) 


0.16390(94) 


0.1633(17) 


0.01 


0.06712(13) 


0.07062(36) 


0.08950(48) 


0.0905(14) 


0.1468(14) 


0.1476(24) 


0. 


0.01503(60) 


0.0250(13) 


0.05746(87) 


0.0589(27) 


0.1330(18) 


0.1341(30) 



13 = 6.10 














0.10 


0.19007( 9) 


0.19200(10) 


0.19945(14) 


0.19986(21) 


0.21820(28) 


0.21745(27) 


0.05 


0.13006(10) 


0.13114(19) 


0.13759(25) 


0.13820(29) 


0.16415(33) 


0.16414(34) 


0.03 


0.09990(10) 


0.10070(21) 


0.10733(20) 


0.10767(34) 


0.14073(39) 


0.14089(53) 


0.02 


0.0816(11) 


0.08239(22) 


0.08985(21) 


0.08957(41) 


0.12870(47) 


0.12837(75) 


0.01 


0.0584(13) 


0.05924(32) 


0.06889(20) 


0.06910(57) 


0.11589(72) 


0.1134(13) 


0. 


0.0112(10) 


0.0135(19) 


0.03684(47) 


0.0360(13) 


0.10381(87) 


0.1012(16) 



curve fitting^^) to the fit of tlie correlator to the form H2-18p . In principle the con- 
strained curve fitting enables a fit of correlators to multipole forms with arbitrary 
number of terms. In practice, however, we find that the mutipole fitting does not 
produce a stable result partially because of the large fluctuations. We therefore flt 
the data to the form (|2-18|) in relatively narrow flt range where excited modes are 
sufficiently reduced according to the observation of the effective mass plot. 

The result for the spectrum is listed in Table [VJ The quoted errors are statistical 
errors evaluated with jackknife method. The meson masses for the four lightest quark 
masses are displayed in Figure [Tj The figure shows that tt and tt, tts and tts, Pq and 
Pq are respectively degenerate. All pionic channels approach to the masses of vr, the 
Goldstone pion channel, as /? increases. This behavior is a signal of flavor symmetry 
restoration toward the continuum limit. The size of flavor symmetry breaking effect 
is estimated with the mass of, say, tt^. Our result of the tt^ mass is the same size as 
the result of Ref.^^-* at /3 = 6.0 (a ~ 2 GeV) on an isotropic lattice. 

To determine the masses at the chiral limit, we fit the lightest three masses in 
each channel to the following forms: 

IT^li^q) = "1-^(0) + bninq (tT, TT, VTs, TTa), (5-1) 

mp{mq) = mp(0) + bpTUq {p^, p^). (5-2) 

The results of the fits are also listed in Table IVI and shown in Figure [7| The above 
forms well represent the data. We also perform a fit of meson mass squared to a 
quadratic form in niq using the four smallest meson masses in each channel. The 




Fig. 7. Meson spectra at P — 5.75 (top-left), 5.95 (top-right), and 6.10 (bottom- left). The curved 
lines represent the results of linear fits using the data at the lightest three quark masses. The 
bottom-right panel shows the meson masses in the chiral limit in physical units. 



results of the latter fits are consistent with those with the former linear fits. 

In the Goldstone pion channel, a small mass remains even in the chiral limit. This 
is considered due to the finite volume effect and the uncertainty in 7^. The masses 
of other pionic channels in the chiral limit decrease as /3 increases. This behavior 
is shown more apparently in the bottom-right panel in Figure [71 which displays the 
meson masses in the chiral limit in physical units as the function of a^. Since the 
staggered action contains 0(0^) systematic uncertainty, the flavor symmetry break- 
ing should disappear linearly in a^. Our result is consistent with this prediction, 
while the finite lattice artifact is sizable even at our largest /?. Rather large lat- 
tice artifact is also expected from the discrepancy between two calibration schemes 
examined in this work. The p meson channel is also consistent with the expected 
behavior. The masses of p in the chiral limit seems to approach the experimental 
value, nip = 770 MeV, within the 0(10%) error of the quenched approximation. 
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To investigate how the uncertainties of affect the masses, we evaluate 

Rn = ^-^H^eI^ (5-3) 

assuming linear dependence of mn in 7^7 in the vicinity of 7^. The subscript H 
specifies the channel. As a general tendency, R increases as quark mass decreases. 
This means that the accurate determination of 7^ becomes increasingly important as 
quark mass decreases. R is in general negative. In pionic channels the absolute values 
of R are 1-2 at /3 = 5.75 and 5.95 and around 1 at /? = 6.10 for ruq = 0.01, while 
at niq = 0.1 it decreases less than half the values at niq = 0.01. By adding the two 
errors in Eq. (|4-3() in quadrature, the error of 7^ is evaluated as 0.11, 0.12, and 0.07 
for (3 = 5.75, 5.95, and 6.10, respectively. At rriq = 0.01 therefore the uncertainty 
in the meson mass due to the uncertainty in 7^ amounts to 20% at P = 5.75 and 
5.95, and 7% aX (3 = 6.10. These uncertainties explain the nonvanishing masses 
of Goldstone pion at the chiral limit. For the vector channels, the values of R are 
about 1.5 times larger than the pionic channels. These quite significant effects call 
for precise calibrations in order to use the anisotropic staggered quark action in 
practical simulations. 

§6. Conclusion 

In this paper, we performed a calibration of the staggered quark action on 
quenched anisotropic lattices with the renormalized anisotropy ^ = 4. As the cali- 
bration procedures, we adopted two schemes which define the fermionic anisotropy 
through the masses in the fine and coarse directions, and through the meson dis- 
persion relation. At values of (3 explored in this work, these two schemes produce 
inconsistent result even in the chiral limit except for our highest (3 {a^ — 2 GeV). 
Although the discrepancy seems to disappear toward the continuum limit, these 
rather large discrepancy should be regarded as a source of significant systematic 
uncertainties. 

In the second part of this paper, we performed a meson spectroscopy with the 
wall source using the values of 7^ determined in the calibration. The result is 
consistent with the result on an isotropic lattice. ^^-^ Toward the continuum limit, 
the flavor symmetry breaking effect seems to disappear. However, the uncertainty 
of 7i? causes severe effect on the meson masses when the quark mass approaches 
the chiral limit. For precision studies using staggered anisotropic quarks, precise 
calibration is indispensable. 

This work displayed how the staggered quark action is realized on anisotropic 
lattices in the quenched approximation. Except for the rather large artifact signaled 
by the discrepancy between the two calibration schemes, no unreasonable result was 
found. However, there are numbers of subject which should be examined still in 
the quenched simulations. It is important to study in more detail how the discrep- 
ancy between the calibration schemes disappear as one approaches the continuum 
limit and to quantify the effect on the observables such as the hadron masses and 
decay constants. The other important subject is to examine how the various chan- 
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nels behave when they are measured in the coarse direction or boosted into a finite 
momentum state. These studies are important to investigate the structure of the 
staggered quark formulation as well as for applications to the calculations for which 
the fine temporal cutoff is crucial. For intermediate range of lattice cutoffs, develop- 
ment of improved staggered quark action is also important to avoid the large lattice 
artifact. The results of this paper provide fundamental information for such further 
studies. 

§7. Toward dynamical simulations 

In this final section, we discuss what we learn from the quenched results in this 
work toward the dynamical simulations. Since in dynamical simulations the gener- 
ation of gauge configurations requires resources, one wants to finish the calibration 
as quickly as possible. In this sense, the mass ratio scheme is economical since the 
statistical error in S,f is smaller than the dispersion relation scheme. However, larger 
quark mass dependence of 7^ is inconvenient to survey the global dependence of 7^ 
on (3 and ruq. In this sense, the dispersion relation scheme has an advantage. There 
is another advantage in the latter: one needs to have a large extent only in the tem- 
poral (fine) direction. The spatial direction can be kept modest size with which the 
momentum modes are not very large. One can also use the extended source field to 
reduce the statistical fluctuations. For these reasons, we realize that the dispersion 
relation scheme is more suitable for the calibration of dynamical lattices. 

It is important to forecast 7^ and to avoid waste of simulation. The mean- 
field estimate is helpful for this purpose. At the quenched level, the mean-field value 
provides a good estimate of the anisotropy parameter of the gauge action. In this 
work, we showed that the tuned anisotropy 7^ approaches to the mean-field value 
from below as P increases. These may help us to select the values of 70 and 7^ to 
be explored. 

Our preliminary result of the calibration of dynamical anisotropic lattices with 
the same combination of actions as in this work indicates that the renormalized 
anisotropy defined through the gauge observable less depends on 'jp, and one through 
the fermionic observable less depends on 7g-^^'* These dependences of and on 
7G and indicate that linear fits may suffice to determine 7^, and 7^ satisfying 
?g(7g'i7f) ~ Ci^(7G'7F) ~ ^- would be able to perform the calibration by 
generating configurations at best at four (or three) sets of (70, 7f)- 
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